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1.  Introduction. 

Of  primary  interest  in  this  project  are  systems  that  are 
modelled  as  beams  or  thin  plates  which  may  undergo  either  small  or 
large  deformations.  Hence,  linear,  so  called  Kirchhoff,  and 
nonlinear,  so  called  von  Karman,  are  amoung  the  model  equations  to 
be  studied.  Linear  models  are  derived  under  assumptions  of  pure 
bending  with  deformations  that  are  small  compared  to  thickness 
which  in  turn  is  assumed  to  be  small.  Larger  deformations  give 
rise  to  in  plane  forces  from  the  stretching  of  the  middle  surface 
and  lead  to  nonlinear  models,  [22,23].  For  beams  linear  models 
with  various  linear  and  nonlinear  damping  terms  are  among  those  of 
interest . 

An  identification  problem  seeks  to  determine  parameters 
within  a  mathematical  model  from  observed  data.  The  central  issue 
is  how  to  utilize  data  to  determine  the  desired  parameters  with 
in  the  context  of  the  model.  Thus,  available  is  a  model  equation 

(1.1)  L(q)u(q)  =  f (q) 

in  which  the  parameters  q  to  be  estimated  belong  to  a  specified 
admissible  set  Qa(j  of  a  Banach  or  Hilbert  space  Q.  The  set  Qa(j 
should  be  physically  meaningful,  such  that  the  mapping  q  •— *  u(q) 
is  well-defined  from  the  state  space  X,  and  such  that  q  *—  u(q) 
is  continuous  with  respect  to  suitable  topologies  on  Q  and  X. 
Available  are  data  z  which  are  in  the  form  of  measurements  on  the 
system.  We  view  these  data  as  belonging  to  an  observation  space  Z 
(Hilbert  space).  One  seeks  to  determine  a  parameter  q  such  that 
Cu(q)  "matches"  z  in  a  suitable  sense.  Here  C  represents  an 
observation  operator  that  maps  X  into  Z. 
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One  approach  to  this  problem,  sometimes  called  the  "regular¬ 
ized"  output  least  squares  method,  cf [ 1-6 , 9-19 , 2 5- 3 5 ] ,  is  formu¬ 
lated  as  a  minimization  problem: 

Find  qg  «  Qa^  such  that 

(1.2)  J(q0)  =  inf imum{ J(q) :  q  s  Qad) 
where 

(1.3)  J(q)  =  IICu(q)  -  z!l|  +  e  ?!q!!^ 

with  p  l  0.  This  formulation  is  given  with  Q  and  Z  Gilbert 
spaces.  Solving  this  problem  often  requires  constraints  on  Qad 
to  provide  (i)  existence  of  a  solution  to  the  state  equation 
(1.1)  and  (ii)  sufficient  compactness  to  obtain  existence  of  a 
solution  to  (1.2)— (1.3) .  Issues  of  interest  for  (1.2)— (1.3) 
include  existence,  regularity  of  the  solution,  approximation, 
stability  with  respect  to  data  or  constraints,  and  uniqueness  or 
identif iabilty. 

If  one  expects  to  use  estimation  techniques  to  identify 
parameters  from  measurements,  then  one  must  first  decide  upon  the 
basic  form  of  the  mathematical  model  (1.1).  Hence,  one  must 
decide  whether  the  mathematical  model  is  appropriate  for  the 
physical  system  being  observed.  The  mathematical  model  embodies 
the  pertinent  principles  a. id  assumptions  from  physics  and  contin¬ 
uum  mechanics  in  addition  to  those  made  on  the  basis  of  geometry, 
etc.  Certainly,  one  must  be  cautious  about  using  models  to  fit 
data  that  do  not  satisfy  the  assumptions  upon  which  the  models 
are  based. 

Having  chosen  the  mathematical  model,  one  must  address  how  to 
use  the  data  within  the  framework  it  imposes.  The  data  for 
systems  of  interest  in  this  work  are  assumed  to  be  obtained  as 
pointwise  measurements  of  deformation,  velocity,  or  acceleration 
or  their  Fourier  transforms.  These  data  typically  should  be 
processed  in  some  manner  and  perhaps  smoothed  before  their  use. 
Also,  data  should  be  in  a  form  consistent  with  the  mathematical 
model  having  the  proper  units,  etc.  The  mathematical  analysis  of 
identification  problems  such  as  (1.1)-(1.3)  guides  in  how  best  to 
formulate  estimation  to  use  these  data.  Of  direct  interest  to 
this  analysis  is  to  determine  how  to  design  and  collect  data  to 
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utilize  desirable  theoretical  properties.  Conversely  if  data  is 
available  from  a  system,  an  important  issue  is  to  determine  how 
much  one  may  reasonably  expect  to  determine  from  the  data. 

In  Section  2  ve  indicate  results  that  ve  have  obtained  on 
the  estimation  of  elastic  coefficients  in  structures  composed  of 
beams.  In  Section  3  ve  report  our  results  on  the  estimation  of 
elastic  parameters  in  nonlinear  static  models  of  thin  plates 
undergoing  large  deformations.  In  Section  4  ve  consider  the 
estimation  of  electrical  conductivity  from  laboratory  measurements 
of  electrical  potential.  This  is  of  interest  since  in  this  study 
ve  have  available  laboratory  measurements  vith  vhich  to  test  our 
algorithms.  Finally,  in  Section  5  ve  indicate  some  results  on  the 
stability  of  sets  of  optimal  estimators  vith  respect  to  data  as 
vieved  as  set-valued  mappings. 


2.  Estimation  in  a  connected  beam  model. 

During  this  project  ve  have  initiated  an  investigation 
concerning  the  estimation  of  coefficients  or  the  design  of  beam 
systems.  In  particular,  ve  have  studied  the  folloving  simple 
problem.  Suppose  tvo  beams,  1  and  2,  of  length  !?j_  and  '?2 , 
respectively,  are  joined  at  right  angles  to  one  another.  The  other 
end  of  beam  1  is  clamped  and  the  remaining  end  of  beam  2  is  free. 
In  response  to  a  force  that  is  perpendicular  to  the  plane 
containing  the  beams  are  deformations  v^_  and  V2  of  beams  1  and  2, 
respectively,  and  a  rotation  9  of  tvist  along  beam  1.  For  the 
static  problem  ve  study  the  folloving  system  of  elliptic  equations 
vith  coupled  boundary  conditions  in  vhich  x  and  y  represent  local 
coordinates  for  beams  1  and  2  respectively. 

(a  wxx>xx  =  fl  in  (0'  V 
( 2 . 1 ) ( i )  -(b  9X)X  =  g  in  (0, 

(c  w2yy ) yy  =  f2  in  {0'  °2 1 

vith  boundary  conditions 

vj^O)  =  vlx(0)  = 


0 
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e(0)  =  o 

at  the  clamped  end, 

(  2 . 1 )  ( 1  i  )  wlxx(  !?1  >  =  0 

(2.1)  (iii)  wx(  )  =  w2(0) 

(2.1)  (  iv)  0(  '\)  =  w2y(0) 

(2.1)  (v)  (b  ex) ( Dj.)  =  (c  w2yy) (0) 

( 2 . 1 )  ( vi )  (a  vlxx)x(  ^ )  =  (c  w2yy ) y ( 0 ) 

at  the  junction,  and 

w2yy(,'2)  =  ( c  w2yy)y(l?2)  =  0 

at  the  free  end.  Equations  (2.1)  are  the  Euler  equations  for  the 
potential  energy  functional 

(2.2)  p  ( v  ,  0,  v2 )  =  JQ1  (a(x)  v2xx(x)  +  b(x)  9x(x)l  dx  + 

+  Tq2  c(y)  v|yy(y)  dy  -  2  j^1  tf^fx)  Vj^x)  +  g(x)  0(x)]  dx  - 

-  2  ,f02  f2(y>  w2(y)  dy- 

The  solution  of  (2.1)  may  be  obtained  by  minimizing  the 
functional  (2.2)  over  the  Hilbert  space 

V  =  {v  =  (v1,e,w2)  «  H2(0,'?1)  x  H1(0,!?1)  x  H2(0,Q2)  : 

vL(0)  =  wlx(0)  =  0,  6(0)  =  0, 
e(  <?1 )  =  v2y(0),  vx(  '^)  =  v2(0)>. 

To  solve  (2.1)  numerically,  the  space  of  basis  functions 
should  satisfy  the  conditions  at  the  junction.  However,  it  is 
much  more  difficult  to  enforce  the  essential  boundary  conditions 
at  the  junction  on  the  basis  elements  that  are  used  to  construct 
the  finite  dimensional  approximating  space  as  one  does  for  example 
with  the  clamped  boundary  conditions.  We  take  two  approaches  to 
this  problem.  The  first  is  to  use  basis  functions  that  are  non- 
conforming  at  the  junction  and  hence  do  not  satisfy  the  essential 
boundary  conditions  at  the  junction.  The  conditions  are  enforced 
as  constraints  on  the  discrete  version  of  the  problem  to  minimize 
the  potential  energy  functional.  We  use  a  penalization  or  an 
augmented  Lagrangian  method  to  impose  the  conditions.  The  second 
approach  involves  directly  changing  variables  by  using  the  junc¬ 
tion  constraints  to  reduce  the  number  of  variables  of  the  problem. 
In  this  method  we  solve  for  certain  variables  in  terms  of  the 
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others  by  using  the  junction  conditions.  In  this  way  the  mini¬ 
mization  problem  for  the  potential  energy  becomes  an  unconstrained 
problem  in  the  new  variables.  In  effect  these  new  variables  are 
coefficients  of  basis  functions  that  satisfy  the  junction  condi- 
t ions . 

We  are  confident  at  this  point  that  our  model  is  reasonable. 
This  is  based  on  experimental  spectral  data  that  we  have  obtained 
with  the  aid  of  Prof.  D.L.  Russell  at  MIPAC.  Specifically,  we 
measured  the  e igenf requenc ies  for  an  aluminum  carpenter's  square 
that  was  clamped  at  one  end.  We  obtained  the  first  six  frequencies 
as 

3.75,  10.62,  39.37,  75,  126.87,  and  198.75  Hz. 

Based  on  our  model  using  a  direct  method  for  the  treatment  of  the 
boundary  condition  at  the  junction,  we  obtained  using  a  cubic 
spline  based  scheme  frequencies  of 

2.28,  14.28,  39.69,  78.45,  129.93,  and  194.82  Hz. 
Considering  experimental  and  discretization  error,  we  feel  that 
these  numbers  compare  favorably.  A  new  set  of  experiments  will 
be  run  again  this  spring  with  improved  clamping  apparatus  to 
obtain  a  new  set  of  data.  In  addition  to  comparing  with  experi¬ 
mental  data,  we  have  also  compared  our  results  with  the  NISA 
engineering  package  for  static  deflection  problems  with  point 
loads.  Predictions  for  the  model  equations  (2.1)  agree  to  within 
5  to  6  decimal  places  to  those  from  the  NISA  package. 

We  have  considered  estimation  problems  for  this  structure.  A 
report  of  results  is  in  the  Proceedings  of  the  1989  IEEE  Confer¬ 
ence  on  Decision  and  Control  [33].  In  that  work  we  formulate  a 
sequence  of  estimation  problems  where  in  each  case  the  underlying 
system  is  one  obtained  from  the  augmented  Lagrangian  functional 
associated  with  the  minimization  of  the  potential  energy 
functional.  Hence,  a  sequence  of  estimation  problems  is  defined 
on  the  approximating  systems  obtained  from  the  augmented 
Lagrangian  and  penalty  approximating  scheme.  A  manuscript 
describing  further  results  for  beams  at  different  angles  is  in 
preparation . 


3.  Estimation  of  parameters  in  plate  models. 
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In  work  funded  under  this  grant  we  studied  the 
estimation  of  both  damping  and  elastic  parameters  in  linear 
equations  that  modeled  small  deformations  of  thin  plates  [25-33]. 
In  this  work  we  obtained  results  concerning  the  properties  and 
regularity  of  solutions,  the  approximation  of  solutions,  ?nd  the 
numerical  treatment  of  sample  problems.  In  addition  we  studied 
the  stability  of  regularized  output-least-squares  optimal  esti¬ 
mators  with  respect  to  data,  second  order  sufficient  conditions, 
and  the  error  of  solutions  in  terms  of  sampling  density  of  the 
data  [34]. 

We  wish  to  indicate  in  more  detail  the  investigations  we  have 
been  conducting  on  model  equations  modelling  thin  plates  but  with 
large  deformations.  We  describe  our  work  on  time  independent 
problems.  Currently,  we  are  working  on  time  dependent  problems. 

Our  approach  currently  is  to  use  weak  solutions.  However,  we  are 
also  considering  classical  solutions  as  in  [21,24]  with 
observation;,  before  blow-up  times.  In  the  derivation  of 
approximate  plate  theories,  linear  mathematical  models  are 
obtained  by  including  bending  terms  only.  The  retention  of  only 
these  terms  embodies  the  assumption  that  deformation  is  in  fact 
small  when  compared  to  the  thickness  of  the  plate  which  in  turn  is 
assumed  to  be  small.  For  larger  deformations  it  is  necessary  to 
include  terms  that  model  the  stretching  of  the  middle  plane  [7,8]. 
These  models  give  rise  to  the  so  called  von  Karman  equations. 

(3.1)  Au  =  s  [  .a  ,  u  ]  +  f 

in  £2 

Bo'  =  - [  u ,  u  ] 

where  the  bracket  term  is  given  by 

['*,  i‘)  =  txx  Vyy  +  +yy  +xx  -  2  +xy  +Xy  • 

The  operators  A  and  B  are  typically  fourth  order  linear  elliptic 
operators.  The  function  u  represents  the  deformation  of  the 
plate  from  an  equilibium  position  and  o  is  commonly  called  the 
Airy  stress  function.  Equations  (3.1)  are  accompanied  with 
boundary  conditions,  for  example,  homogenous  Dirichlet  boundary 
conditions 

u  =  iTn  =  0  =  dn  =  0  °n  r. 

We  have  recently  considered  the  estimation  of  the  parameter  a  = 
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a(x,y)  in  the  specific  cases  of 

A  =  -(a  f) 

and 

A  •+'  =  v  ( -  (a  -f )  )  + 


with 


f  ( 1  v )  (  ( a  t'xx  )  xx  +  ^afyy^yy  +  2  (a^Xy^xy^ 


B  '+'  = 


Such  operators  arise  with  variable  coefficients  for  thin  plates 
of  uniform  thickness  but  with  variable  Young's  modulus.  It  turns 
out  that  solutions  u  of  (3.1)  are  stationary  points  of  a  quartic 
functional.  Moreover,  it  is  possible  that  equation  (3.1)  has  more 
than  one  solution  (7,8.203,  and  therefore  the  parameter  to  state 
mapping  may  not  be  well-defined.  Even  so  we  have  demonstrated 
existence  of  an  optimal  parameter  for  the  following  regularized 
output-least  squares  problem. 

(3.2)  Find  aQ  •s  Qad  such  that  J(ag)  =  inf{J(a):  a  s  Qa(j} 


where 


J(a) 


!!u  ( a  )  —  z  !!l  i  (  Q)  *  £  1 1  a  1 1  h  £  {  q  j 


(3.3)  Qad  =  (a  •=  H2(Q)  :  a  l  S  >  0}  . 

If  there  may  exist  multiple  solutions  to  the  equation  (3.1),  the 
fit-to-data  functional  above  is  not  well-defined.  However,  it  can 
be  shown  that  the  solution  set  U(a)  of  equation  (3.1)  is  closed 
and  compact  in  the  weak  topology  of  Hg(£2).  Since  the  embedding  of 

H  q  ( C2 )  into  L^(Q)  is  compact  when  Q  is  an  open  domain  in  P?  with  a 
Lipschitz  boundary,  we  denote  by 


!lu(a)  -  zllL2(H)  =  min{hU  -  zi!L2(,-.):  u  s  U(a)  for  a  s  Qad3. 

The  minimimum  exists  by  the  continuity  of  the  L^-norm  with 
respect  to  the  weak  topology  on  H§(£).  It  can  be  shown  that 

there  exists  a  solution  to  the  problem  (3.2)— (3.3) . 

Let  us  set 

il-t'ilo  =  :i'f!!Li  and  ll'fll  *  II^’+'IIq. 

It  is  well-known  that  thus  defined  is  a  norm  on  Hq(S2).  More¬ 

over,  there  is  a  constant  k  such  that  for  any  ~  Hq(Q), 

’ ! 'F 1 1 o  i  k  .  In  addition  there  is  a  constant  K  such  that  for  u. 


:'v!i . 


v,  and  v  in  Hg(Q) 

I  jo  [U,  v]  w  dx  !  <  K  Hull  Hv'l 

Set  p  =  (k  II  f  !l  o )  /  3  and  r  =  e  K2/  ’. 

It  is  natural  to  consider  linear  approximations  to  (3.1). 

By  imposing  conditions  on  e,  v,  and  f  we  may  obtain  by  the 
contraction  mapping  principle  appr ox ima t i ng  linear  systems  for 
which  convergence  may  be  established.  Based  on  this  approach  we 
obtain . 

Theorem .  Let  M  be  chosen  such  that  ^  <  p  <  M  and  assume  r 

satisfies  r  <  (M  -  ,p>/(M').  Then  for  any  a  e  Qad  there  is  a 
unique  solution  u(a)  of  (3.1)  with  !!u(a)l!  *  M. 

In  (311  we  consider  the  following  approximating  linear 
systems . 

(3.4)  U_^  =  i?g  =  0 

for  i  =  0,  1,  .  .  . 


-( a  ^  u^  =  e  +  f  in  Q 

ui  =  Jni  =  0  on  r 

(3.5) 

“2  0 i + l  *  ~  ui^  in  - 

ji  =  0  on  r- 

Theorem .  With  M,  -,  and  r  defined  as  above,  the  sequence  u ^ 

generated  from  (3.4)— (3.5)  converges  uniformly  to  u(a)  in  H^(Q) 
for  a  •=  Qacj  where  llu(a)ll  i  M  and 

l!u  ±  -  u  ( a  )  I!  <  M  (  3  r“M)i/(l-3riM)  . 

The  assumptions  that  2M/3  <  F  and  r  <  (M  -  ?)/M7  imply  that 
3r  "M  <  1.  Hence,  if  we  consider  the  mapping  a  *  u(a)  of  Qa(j 

into  the  ball  (u  s  H^(Q):  Hull  i  M>,  then  it  is  well-defined.  We 
may  demonstrate  differentiability, of  solutions  of  (3.1)  with 
respect  to  a  under  the  assumption  that  3r~M  <  1.  Therefore, 
regularity  results  for  optimal  estimators  may  be  obtained  under 
the  above  conditions. 

Theorem .  Let  3r “M  <1.  If  a  is  a  solution  of  (3.2)— (3.3),  then  a 
H2+n(Q)  for  n  e  (0,1). 

These  results  enable  us  to  provide  an  approximation  theory  for 
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these  problems. 

We  base  an  estimation  algorithm  on  the  above  approximating 
problems  in  which  we  estimate  the  parameter  a  in  each  linear 
approximating  system  and  use  the  computed  state  for  the  update  in 
the  next  linear  problem.  Details  and  the  results  of  numerical 
experiments  are  reported  in  [31]. 

The  above  formulation  imposes  rather  strong  conditions  on 
f,  z  and  :•  in  order  to  apply  the  contraction  mapping  princple. 
These  conditions  also  imply  differentiability  of  the  mapping  a  1 — 
u(a).  We  note  however  that  existence  of  optimal  estimators  does 
not  depend  upon  the  these  conditions.  We  consider  the  following 
weak  model  error  formulation.  Define  w  :  H2(£i)  x  V  *  V  as 
follows 

(3.6)  -2  w  =  -(a  o.  u)  -  z  [a,  u]  -  f 

-2  a  =  -  [  u ,  u  ]  i  n  '.J 

with  boundary  conditions 

w  =  dn  =  0  on  r 

It  is  easy  to  see  that 

(3.7)  w (a, u )  =  0 

if  and  only  if  equation  (3.7)  is  satisfied.  We  can  show  that  the 

equation  (3.7)  implicitly  defines  a  function  a  '  u(a)  if 

_1 

(3.8)  1  -  K  ( 4 )  “  ilf  I'h_2(Q)  >  0 

and  this  function  is  differentiable.  Thus,  we  may  weaken  condi¬ 
tions  to  obtain  local  existence  and  local  differentiability  of  the 
parameter -to -state  mapping. 

The  optimization  problem  in  this  formulation  is  given  as 
Find  ( a q , u q  )  such  that 

J(a0,UQ)  =  inf{J(a,u):  a  s  Qad  and  w(a,u)  =  0} 

where 

J(a,u)  =  !Iu  -  zll2  +  ^  ,!a'!^i. 

Again  under  the  condition  (3.8),  we  may  obtain  the  existence  of 
Lagrange  multipliers  associated  with  the  constraint  (3.3).  This 
Lagrange  multiplier  may  be  estimated  in  terms  of  llug-zily.  Thus 
we  may  find  conditions  for  the  choice  of  ?■  and  the  positive 
definiteness  of  the  second  derivative  of  the  Lagrangiar.  func¬ 
tional.  With  these  estimates  we  may  obtain  local  uniqueness  and 
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stability  results  and  show  that  the  augmented  Lagrangian  method 
converges  [32].  These  results  have  been  submitted  for  publi¬ 
cation.  They  have  also  been  presented  at  the  Midwest  Differ¬ 
ential  Equations  Conference  in  November  1989. 

Without  further  assumptions  such  as  those  in  (3.7)  ie  may 
penalize  the  constraint  (3.6).  This  general  formulation  is  given 
by  means  of  the  minimization  problem 

Minimize  ! I u  -  zlly  +  =-  +  7  !lw(a,u)  (ijjji 

We  have  analyzed  this  problem  and  have  presented  initial  results 
in  a  colloquium  at  the  Division  of  Applied  Mathematics  at  Brown 
University  last  spring. 

We  include  some  preliminary  results  of  numerical  tests  using 
both  the  penalty  method  and  the  augmented  Lagrangian  for  this 
problem.  In  this  case  we  consider  the  equation 

-2  v(a,u)  =  -(  a  -u)  +  £  [B(u,u)/u]  -  f 
with  w  =  w(a,u)  s  Hg(Q)  and 

-2  B(u,u)  =  [  u,  u ] . 

We  consider  the  minimization  problem: 

Minimize  L(a,u;')  subject  to  u  ■=  V  and  a  •=  H2 

for 

L(a,u;  )  =  7-  llu  -  z!l27  +  4  Hall2?  +  (x,w(a,u))v  +  £  II  w  ( a ,  u )  II 2 

4  L4  c  v  4  y 

where  V  =  Hg.  We  generate  f  by  determining  f  from 

f  =  ^  (  aT  -uT)  +  s  [B(uT,uT) ,uT] 

where  uT  =  256  x2  y2  (1-x)2  (1-y)2.  We  used  discretizations  with 
8  subintervals  in  both  x  and  y  directions  for  approximating  u  and 

3  subintervals  in  both  x  and  y  directions  for  approximating  a. 

Accordingly,  the  mesh  for  the  example  that  we  give  is  very  coarse. 
Approximating  functions  were  tensor  products  of  cubic  B-splines 
adjusted  for  boundary  conditions.  We  used  the  conjugate  gradient 
method  for  the  optimization  steps.  Our  initial  guesses  were  off 
by  relative  L2  errors  of  87%  and  45%  for  u  and  a,  respectively. 
Below  we  give  values  of  the  exact  and  computed  functions  for  the 
state  u  and  the  coefficient  a.  As  an  example  consider  a  problem 
with  a  discontinuous  elliptic  coefficient 
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fl.5,  (x,y)  such  that  0.25*y*0.75 
aT  =  |1.0,  otherwise 

with  initial  relative  L2  errors  of  79.4%  and  231%,  for  u  and  a, 
respectively.  After  19  iterations  we  obtained  the  following  with 
the  augmented  Lagrangian  :«ethod. 


■  (0,0) 

(  0 . 2 . 0 . 2  ) 

O 

* 

rr 

O 

(  0 . 6 . 0 . 6  ) 

(  0 . 8 . 0 . 8  ) 

Uip  |  0 

0 . 168 

0.849 

0.849 

0.168 

ucalc  1  0 _ 

0.111 

0 .989 

0.849 

0.260 

For  the  coefficient  a  we 

have 

O 

* 

o 

(0.2. 0.2) 

(0.4.0. 4) 

(0.6. 0.6) 

(0.8. 0.8) 

3rj»  j  X  •  0  0 

1.00 

1.50 

1.50 

1.00 

acalc  !  0-649 

0.948 

1 .67 

1.67 

0.925 

II  w  ( a ,  u  )  llL  2  =  0.0884. 

For  the  penalty 

method  we 

have  for  the 

state 

o 

o 

(0.2. 0.2) 

(0.4. 0.4) 

(0.6. 0.6) 

(0.8.0. 8) 

Ut  j  0 

0.168 

0.849 

0.849 

0.168 

ucalc  i  0 

0.171 

0 .844 

0.794 

0 . 0893 

For  the  coefficient  we  obtain  the  following. 

o 

k 

o 

(0.2. 0.2) 

(0.4.0. 4) 

(0.6. 0.6) 

(0.8. 0.8) 

o 

o 

H 

1.00 

1 . 50 

1 . 50 

1 . 50 

acalc  i  0-533 

0.603 

1.43 

1.51 

0.619 

Ilw  ( a ,  u  )  il  L  2  =  0.140. 


4.  Estimation  with  experimental  data. 

In  the  previous  proposal  we  expressed  one  goal  as  that  of 
considering  estimation  problems  using  data  that  has  been  obtained 
experimentally.  We  feel  that  this  is  the  essential  next  step  in 
producing  useful  estimation  algorithms.  We  report  our  efforts 
here  for  an  experiment  to  obtain  electrical  conductivity  from 
potential  measurements  in  a  circular  reservoir.  The  situation  is 
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the  same  as  that  of  determining  elasticity  coefficients  for  a 
membrane  from  displacement  measurements.  It  is  shown  in  [8]  that 
such  a  membrane  equation  may  be  obtained  as  a  limiting  case  of  von 
Karman  equations. 

The  experimental  setup  was  as  follows.  A  cylindrical  tank 
with  circular  cross-section  Q  is  filled  with  sand  in  which  there 
is  a  salt  water  solution.  There  are  eleven  locations  in  7:  at 
which  a  source  electrode  may  be  located  or  potential  may  be 
measured.  There  is  a  ground  located  at  an  additional  position. 

The  walls  of  the  tank  are  nonconducting.  The  electrodes  that  are 
used  are  such  that  the  current  is  focused  so  as  to  be  two  dimen¬ 
sional.  This  allows  us  to  consider  a  two  dimensional  problem.  At 
a  certain  point  in  time,  a  volume  of  salt  water  is  extracted  and 
replaced  by  injecting  fresh  water.  Consequentially,  the  electi- 
cal  conductivity  is  a  spacially  dependent  function.  An  electrical 
source  of  current  is  placed  at  a  position  and  measurements  are 
then  made  at  the  remaining  locations.  These  measurements  are 
repeated  with  the  source  at  various  locations.  The  diffusion  of 
the  salt  into  the  fresh  region  is  sufficiently  slow  as  compared 
with  the  time  it  takes  to  make  the  measurements  that  the  problem 
may  be  considered  time  independent. 

The  governing  model  is  given  by  a  two  dimensional  elliptic 
Neumann  problem 

--’(a”u)=fin2 
0  on  32. 

where  u  is  the  electrical  potential  function,  the  coefficient  a  is 
the  electric  conductivity,  and  f  is  the  source  term.  Here  for  the 
purpose  of  mathematical  formulation  we  consider  f  to  be  given  by 

f  =  T  Is 

where  !s  is  a  function  of  small  support  (in  a  ball  centered  at 
the  source  location  xs)  such  that 


,U  5s(x)  dx  =  1 


Measurements  of  the  electrical  potential  are  made  at  Nobs 
locations  over  the  domain  71.  Thus  the  observation  operator  C 


N 


takes  u  into  P 


obs 


is  given  by 
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C  u(a)  =  coH<u(a)/jafj>:i=l/ .  .  .  #Nobs> 
where  0^  are  continuous  linear  functionals  on  H^CCiHor  possibly 
H2(f2)  depending  on  the  regularity).  For  example,  we  can  take 


where  is  a  disk  centered  at  the  observation  location  x ^  to  be 
an  observation  functional.  The  measurements  are  given  as  an 
N0bs-vect:or  z  real  numbers. 

To  treat  the  equation,  we  consider  the  weak  model  error 
w  =  w(a,u)  defined  as  the  solution  of 

-£.v  +  v  =  -~*  (a  "  a)  -  f  in  5 

with 

boundary  condition 

fj£  =  0  on 

The  estimation  problem  is  now  given  in  terms  of  the  optimi¬ 
zation  problem 

Find  ( a  q , u  q  )  s  Q  x  H1(Q)  such  that 

J(a0,ug)  =  inf{J(a,u):  a  «  Q  and  u  •=  H^- ( '.!)  satisfy  w(a,u)  =  0} 
where 

J(a,u)  =  v  ICu  -  zl2  +  p:  lla IIq , 

Q  =  H2(Q),  and  !  •  !  is  the  Euclidean  norm  in  H?N°bs. 

To  solve  this  problem  numerically,  we  approximate  the  disk 
centered  at  the  origin  with  quadrilaterals  obtained  by  partition¬ 
ing  the  y-direction  with  an  odd  number  of  Ny  equally  spaced 
horizontal  lines  where  the  (Ny-  l)/2  +  1  is  line  y  =  0.  Each 
horizontal  segment  we  partitioned  into  Nx  sub interva Is .  Connec¬ 
ting  the  corresponding  points  forms  the  system  of  finite  elements. 
This  discretization  may  be  mapped  to  a  square  region  that  is 
partitioned  similarly  by  defining  transformations  that  map  each 
subquadrilateral  to  the  corresponding  subrectangle.  We  use  a 
system  of  basis  functions  that  are  bicubic  functions  made  up  of 
tensor  products  of  cubic  B-splines  as  basis  functions  for  the 
square.  The  composition  of  the  transformations  with  these  func¬ 
tions  then  determines  basis  functions  on  the  disk.  We  used  64 
basis  functions  to  approximate  the  state  and  49  to  approximate  the 
coefficient  a. 
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The  estimation 

problem  is  solved  by  means  of  the  augmented 

Lagrangian  method. 

Thus  we  form  the  functional 

F  ( a , 

u;  x,K ) 

2 

=  J(a,u)  +  <  a  ,  w  ( a ,  u  )  >  +  K  ilv(a,u)!lHi 

and  use  the  conjugate  gradient  method  for  minimization  at  each 

step.  We  use 

v  =  1 

.0,  r  =  10“  ,  and  K  =  10.0.  The  coefficient 

is  calculated 

for  two  different  sets  of  data  to  compare  the 

estimated  a's 

.  We 

obtain  the  following. 

CASE  1. 

X 

y 

values  of  a  (est) 

-53.2 

-109.8 

2 .87 

-17.7 

-109.8 

3.09 

17.7 

-109.8 

2.99 

-116.4 

-36.6 

2.83 

-38.8 

-36.6 

3 . 04 

38.8 

-36.6 

2.82 

-116.4 

36.6 

2.48 

-38.8 

36.6 

2.80 

38.8 

36.6 

2 .75 

The  relative 

is  4.9%  while 

&  t  y  rw 

^obs 

.Z  (u(xi,yi )  -  )l 

C  L  L  KJ  L. 

l!wllHi 

L  UI  Li  04I  L  \  ) 

Nobs  7 

21  Z; 

i  =  1  l 

=0.071.  In  the  second  case  we  have  the 

following 

CASE  2. 

X 

y 

values  of  a  (est) 

-53.2 

-109 . 8 

2.91 

-17.7 

-109 .  a 

3.11 

17 . 7 

-109.8 

3 . 00 

-116.4 

-36.6 

2.87 

-38.8 

-36.6 

3 .05 

38.8 

-36.6 

2 .83 

-116 . 4 

36.6 

2.51 

-38.8 

36.6 

2.82 

38.8 

36.6 

2 . 78 

with  relative  error  for  u  =  4.4%  while  l!vllHi  =0.038.  Thus,  we 
are  able  to  estimate  the  coefficient  a  so  that,  for  different 
sets  of  data,  the  values  of  the  estimated  function  agree  to 
within  the  accuracy  of  the  data.  We  feel  that  this  is  as  much  as 
one  can  reasonably  expect  and  is  a  successful  application  of  our 
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methods.  This  work  will  appear  as  an  invited  paper  in  a  special 
issue  on  parameter  identification  in  the  Journal  on  Advances  in 
Water  Resouces  [35]. 

5.  Stability  properties  with  respect  to  data. 

We  have  studied  stability  properties  of  optimal  estimators 
with  respect  to  data.  Except  for  results  in  [13]  in  which  there  is 
strong  regularization  and  it  is  assumed  that  the  data  is  in  the 
attainable  set  solutions  of  output  least  squares  estimation 
problems  are  not  unique.  In  [9,31,34]  conditions  are  given  for 
regularization  under  which  there  is  a  certain  degree  of  stability. 
These  results  do  not  require  that  the  data  belong  to  the  attain¬ 
able  set  and  do  not  provide  for  uniqueness.  We  have  thus  begun  a 
study  of  the  continuity  properties  of  the  set  of  optimal  esti¬ 
mators  as  a  set-valued  function.  To  be  more  specific,  we  study  the 
following  sample  problem.  Let  Q  be  a  bounded  open  domain  in  IR11 
with  a  Lipschitz  boundary  n.  We  consider  the  following  sample 
boundary  value  problem 

(5.1)  -'7'(a'?u)=finQ 

u  =  0  on  r 

with  f  «  W~l' 2  ( fj)  .  The  coefficient  a  belongs  to  Wk'2(Q)  for  k  =  l,2 
and  satisfies 

(5.2)  0  <  ,»u  I  a  I  a.e.  in  2 
For  a  fit-to-data  functional  we  use 

(5.3)  J(a;z)  =  llu(a)  —  z  II 2  2  +  ?  Mai  Ij2 

where  we  denote  the  dependence  of  u  on  a  by  u(a)  and  the 
dependence  of  J  on  z  by  J(*;z).  The  regularizing  parameter  P  is 
required  to  be  positive.  Finally,  I • is  the  seminorm 

(5.4)  I  !  'fl  \l  =  i"_  A  ,5,  „  ID**!2  dx. 

K  J  Q  0  <  !■:■:!  i  k 

The  estimation  problem  is  given  as 

(5.5)  Find  aQ(z)  s  Qad 

such  that  J(ag(z);z)  =  [J(a;z):  a  =  Qad) 


where 


Qad  =  [a  -s  Wk  • 2 :  a  satisfies  (5.2)). 
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There  exists  a  solution  to  (5.5)  and  we  denote  the  set  of 
solutions  to  (5.5)  by  Q(z)  and  study  the  continuity  properties  of 
the  set-valued  mapping  z  *— Q(z).  We  show  that  in  general  this 
mapping  is  upper  semicontinuous  with  the  weak  topology  or  and 

the  strong  lP  topology  on  the  domain  space  as  a  set-valu  •'  napping 
in  the  following  sense: 

Given  a  neighborhood  N  of  Q(Zq)  there  exists  a  neighborhood  M 
of  Zq  such  that  Q(M)  £  N. 

By  introducing  stronger  regularization,  for  example  with  n=2 
taking  k=2  or  as  in  the  cases  in  [9,31)  or  taking  Qad  to  be  in  a 
finite  dimensional  space  one  obtains  the  stability  results  of 
(9,31).  These  results  amount  to  weak  lower  semicontinuity  of  the 
mapping  z  H?-  Q(z)  in  the  following  sense: 

Given  an  element  ag  s  Q(zg)  and  a  neigh-  borhood  N  of  ag 
there  exists  a  neighborhood  M  of  Zg  such  that  if  z  s  M  then  Q(z)  ;_ 
N  7*  0.  With  the  two  results  together  we  obtain  that  the  mapping 
z  ■—*  Q ( z )  is  continuous  as  a  set-valued  mapping  with  respect  to 
the  Hausdorff  metric  on  the  collection  of  closed  sets  in  the  weak 
topology  on  bounded  sets  of  Qa(j . 
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